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We show that the quantum relaxation process in a classically chaotic open dynamical system is 
characterized by a quantum relaxation time scale t q . This scale is much shorter than the Heisenberg 
time and much larger than the Ehrenfest time: t q oc g a where g is the conductance of the system 
and the exponent a is close to 1/2. As a result, quantum and classical decay probabilities remain 
close up to values P ~ exp (—•^/g) similarly to the case of open disordered systems. 

PACS numbers: 05.45. +b, 03.65.Sq 
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Recently there has been a considerable interest in the 
statistical properties of the poles of the S matrix in meso- 
scopic quantum dots coupled to conducting metallic leads 
The statistical properties of these poles determine 
the effective life time of particles inside the dot and there- 
fore are directly related to the conductance fluctuations 
and current relaxation inside the dot. In fact, the prob- 
lem of current relaxation in diffusive mesoscopic samples 
connected to leads has been addressed long time ago [|| . 
Recently the interest to this problem was renewed and 
new effective methods based on the supersymmetry ap- 
proach have been developed to study the problem in more 
detail Q. For quasi one-dimensional metallic samples 
the results of HQ] predict that the current in the sam- 
ple, being proportional to the probability P(t) to stay 
inside the sample, will decay, up to a very long time, in 
an exponential way according to the classical solution of 
diffusive equation which describes the electron dynamics 
in disordered metallic samples: P(t) ~ exp (— t/t c ). Here 
t c ~ = N 2 /D is the diffusion time for a system of size 
N with diffusion coefficient D. 

According to the strong deviation of quantum 

probability P q from its classical value P takes place only 
for t > tii where the quantum probability decays as 
P q (t) ~ exp (-gin 2 (£/*#))■ Here, t H = 1/A(h = 1) 
is the Heisenberg time, A is the level spacing inside the 
sample and g — tn/t c — E c /A is the conductance of 
the sample with Thouless energy E c — l/t c . At time 
tji, In Pg(tjj)/ In P(tu) ~ 2. As it was pointed out re- 
cently Q, less strong deviations (ln(P q (t q ) / P(t q )) ~ 1) 
should take place at a shorter time t q ~ t cy fg due to 
weak localization corrections according to equations ob- 
tained in Q] . Up to now these theoretical predictions for 
open systems have not been checked neither by numer- 
ical computations nor by laboratory experiments. Also 
the above results are based on an ensemble averaging 
over disorder and their validity for a quantum dynamical 
system which has one fixed classical limit is not evident. 
The investigation of this problem is also interesting from 



the viewpoint of semiclassical correspondence in a regime 
with exponentially fast spreading of narrow wavepackets 
due to which the Ehrenfest time scale |(| is very short: 
tE ~ hxN/A where A is the Liapunov exponent. 

In this paper we study the quantum relaxation process 
in a dynamical model of quantum chaos where diffusion 
is caused by the underlying classical chaotic dynamics. 
This model, introduced in R|, describes a kicked rotator 
with absorbing boundary conditions (when the momen- 
tum is larger than some critical value). This open sys- 
tem can be considered as a model of light trapped in a 
small liquid droplet with a deformed boundary in which 
the rays, with orbital momentum less than some criti- 
cal value, escape from the droplet because the refraction 
angle exceeds the critical value || . 

Contrary to the standard kicked rotator model || in 
which the matrix of the evolution operator is unitary, the 
absorption breaks the unitarity of the evolution matrix 
so that all eigenvalues move inside the unit circle. In 
other words, each eigenvalue can be written in the form 
A = e~" = exp (— iE — T/2) where T characterizes the 
decay rate of the eigenstate. In this way absorption corre- 
sponds to ideal leads without reflections back to the sam- 
ple. A similar approach, in which coupling to continuum 
was studied on the basis of non Hermitian Hamiltonians, 
has been developed and widely used by Weidenmuller et 
al. (see for example ||). 

In our model the quantum evolution of the wavefunc- 
tion is described by the following quantum map: 



iP = Ut/j = Ve 
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where V is a projection operator over quantum states n 
in the interval (— N/2, N/2). Here, the commutator is 
[n, 9] = —i and the classical limit corresponds to k — > oo, 
T — > while the classical chaos parameter K = kT re- 
mains constant. In the classical limit the dynamics is 
described by the Chirikov standard map km: 
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For the classical computations, in analogy with the 
quantum model, all classical trajectories escaped from 
the interval (—N/2, N/2) are absorbed and never return 
back. Due to this absorption, in the regime of strong 
chaos (K ^S> 1) with one chaotic component (no is- 
lands of stability), the classical probability to stay inside 
the interval (—N/2, N/2) decays exponentially with time: 
P(t) ~ cxp(— 7 C £). The time scale t c — l/j c ~ to is de- 
termined by the diffusion time to required to reach the 
absorbing boundary from the center. Since the diffusion 
rate is D =< An 2 > /At ~ k 2 /2 then 7c = E c ~ l/t D = 
k 2 /N 2 . In order to study the quantum relaxation we 
fixed the classical chaos parameter K = 7 and the ratio 
N/k = 4. In this way the diffusion time to is constant 
when TV — ► oo and this allows to investigate the semi- 
classical behavior. Moreover to >> 1 which justifies the 
diffusive approximation. We note also that the system 
(0) with -N/2 < n < N/2 and V = 1 coupled to open 
leads (T = for |rc| > N/2) had been studied in jio). 
The results obtained there showed that this model has 
universal conductance fluctuations |llf| and other prop- 
erties very similar to mesoscopic metallic samples. 

The numerical solution of the classical problem was ob- 
tained by iterating map (||) for M = 9T0 9 different initial 
conditions homogeneously distributed on the line n = 0. 
The results demonstrate a clear exponential decay P(t) = 
exp(- 7c i - b) with 7c = 0.101882(1), b = 0.17774(5) (see 
Fig.l). This exponential decay shows that for K = 7 
the phase space is completely chaotic without any island 
of stability. Even with such a high number of orbits, 
the classical computations allow to obtain directly the 
probability P(t) with 10 % accuracy only up to the level 
P w 5 • 10~ 8 (i w 165). This limitation is due to statisti- 
cal errors appearing for finite number of trajectories. In 
spite of this, the decay rate 7c can be found with very 
high precision which allows to extrapolate the probability 
behaviour to larger times. 

For the quantum evolution we choose the correspond- 
ing initial condition in which only the level n = is popu- 
lated and we studied numerically the quantum dynamics 
(|l|) for different N. We have found that the quantum 
probability P q (t) follows the classical one up to a time 
t q after which it starts to decay at a slower rate (Fig.l). 
We determined the quantum relaxation time t q by the 
condition \n[P q (t q )/ P(t q )] — 0.1 which corresponds to 
10% deviation. The comparison of quantum and clas- 
sical probabilities is shown in Fig.l. The values of t q , 
obtained in this way, strongly fluctuate with changing 
the system size N as is typical for mesoscopic systems 
(Fig. 2). These fluctuations are satisfactory described by 
a log normal distribution (insert in Fig. 2), but a more 
detailed analysis is required to determine precisely their 
statistical properties. 

To suppress the fluctuations, we average hi P q (t) over 
different system sizes by changing N in a small interval 
SN around a given N. Typically we averaged up to 500 



different N values. This allows to determine the averaged 
ratio R(t) — (\n(P q (t)/P(t))) of quantum to classical 
probability. For all these N values the classical dynamics 
is exactly the same since we kept K — 7 and N/k — 4. 
Then the quantum relaxation time t q (averaged) at 10% 
level, is determined by the condition R(t q ) = 0.1. The 
dependence of R(t) on time, for different N , is shown in 
Fig. 3. It is clearly seen that t q grows as we approach 
the semiclassical limit N — > oo. The rescaling of data 
for R(t) as a function of t/t q shows a satisfactory global 
scaling behaviour of quantum probability (see insert in 
Fig. 3). The dependence of t q on N is shown in Fig. 4 for 
the semiclassical regime 500 < N < 130001. This regime 
corresponds to a variation of conductance g = N/t c in 
the interval 50 < g < 13000. The fit of numerical data 
gives a power law dependence t q ss 0.19 • t c ■ N a with 
a w 0.41. This power remains the same for 5% deviation 
level (Fig. 4). 

Here we propose a qualitative explanation of a value 
based on the fact that in the open system the physics is 
affected not by the level spacing A but by the distribu- 
tion of poles of the scattering matrix S which describes 
the coupling to the leads. These poles are located in 
the complex energy plane and their imaginary parts de- 
termine the decay probability of eigenmodes inside the 
sample. For our model, the poles are simply given by 
the evolution operator ([!]). The eigenvalues of U are dis- 
tributed in a narrow ring of width E c inside the unitary 
circle |Q. This is typical for diffusive samples coupled to 
strongly absorbing leads. As a result, N complex eigen- 
values are homogeneously distributed in a ring of total 
area Ak, E c and the distance between them, in the com- 
plex plane, is 5 w \J E c /N . In the classical limit this 
spacing goes to zero and one obtains a continuous den- 
sity of poles. However, for finite N , the separation of 
poles is finite and can be resolved after a time t q ~ 1/5. 
According to this argument, which is independent of the 
symmetry and dimensionality of the problem, the devi- 
ation between quantum and classical probabilities will 
take place at 

t q w 0.38VWV = 0.38t c% /5 (3) 

where the numerical coefficient has been extracted from 
Fig. 4. The theoretical dependence (ph, which corre- 
sponds to a = 1/2, is different but close to the numerical 
value a = 0.41. We attribute this difference to a not suf- 
ficiently large value of y/g. Indeed, neglecting the values 
with yfg < 15, we obtain a = 0.44 for 10 % deviation 
and a — 0.45 for 5% deviation; these values are closer to 
the theoretical prediction a = 0.5. 

The scale t q can be also explained in a more standard 
way based on weak localization corrections |^,^|. In- 
deed, the quantum interference gives a decrease of the 
diffusion rate l/t c oc D — > D(l — at/tji) where a is 
some constant (diffusion stops at time tu)- As a result 
\n(P q /P) Ri at 2 q /(t c t H ) - 1 and one gets (3). 
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For very large times, the decay of P q (t) is determined 
by the eigenvalue e = E — iT/2 with minimal T — T m i n . 
This asymptotic behaviour should start from some time 
scale t f , which can be estimated in a following way. Sim- 
ilarly to the results obtained for complex matrices Jl2| , 
the eigenvalues e should be distributed in a region with 
smooth boundary in the complex plane with approxi- 
mately constant density of points for r ~ y c . Typically 
this boundary is parabolic near the extremal Tq close to 
r m i„. Due to this, the relative number of eigenvalues dW 
in the interval dT is given by dW/dT ~ y/(T - T )/jc /2 - 
The total probability in the interval <5r = T — Tq is 
W ~ (i5r/7 c ) 3 / 2 . The distance between the two lowest 
values of T can be estimated from the condition W ~ 1/N 
which gives ST ~ y c N~ 2 / 3 . Then tj « QAt c N 2 ^ which 
is much larger than t q . The numerical factor here was de- 
termined from the two lowest values of r for N = 5001: 
r TOi „ = Ti = 0.065309 and nearest T 2 = 0.066203. 
These values were obtained by direct diagonalization of 
the matrix U. The rate of the asymptotic decay of P q (i) 
for N = 5001 (Fig.l) agrees, up to 6 digits accuracy, 
with the above value of T m i n . We note that a typical 
size of fluctuations for poles is 6 and so we expect that 
7c 1 rnin ~ 5 ~ l/t q . Since in our model tj < tn 
and the classical limit is fixed we do not see the quan- 
tum deviations discussed in for t > tg- At large 
g one should average over exponentially large number of 
realizations to observe them. 



than the relative size of quantum cell 1/N. Here, the 
error bar gives the average fluctuation of P q (t q — 254) 
obtained for 77 values of N. 

Also it is interesting to note that the Ehrenfest time 
scale is much smaller than the quantum relaxation 
time t q : tE/t q ~ hiN/yfN « 1. For example for 
N = 1.3 • 10 5 we have t q = 254 while t E = In TV/A « 9.4 
(A m In if/2 « 1.25). This shows that the agree- 
ment between quantum and classical relaxation contin- 
ues for a time scale which is much larger than the time 
of wave packet spreading. However, for t > ts there is 
no exponential instability in the quantum motion P,H . 
As a result, correlation functions of the type C(r) = 
(sin#(t) sin#(£ + r)) which, in the regime of strong chaos, 
decay exponentially in the classical case (In |C| ~ —At), 
in the quantum case decay only during the Ehrenfest time 
scale up to In |C| ~ — In N (ig <C r -C t c ). This is simi- 
lar to what happens in closed (unitary) systems such as 
the kicked rotator jL4j. This example shows that expo- 
nential relaxation is not necessary related to exponential 
local instability and positive Kolmogorov-Sinai entropy. 

After the submission of this paper the scale t q has been 
obtained on the basis of random matrix theory and su- 
persymmetry for kicked rotator with random phases [ jl5| . 
The related results for pre-localized states in closed sys- 
tems were discussed in |L6[. 



The largest value of t q we have numerically obtained 
(at N — 1.3 • 10 5 ) is t q — 254 which corresponds to a 
probability P q (t q ) ~ 5 ■ 10~ 12 . On the other hand, the 
classical simulation with M = 9 ■ 10 9 orbits allows to 
directly compute the classical P(t) with 10% accuracy 
only up to t « 165 which corresponds to a probability 
P w 5 ■ 10 -8 . In order to reach the above level of quan- 
tum accuracy one should iterate M w 10 14 orbits, which 
is already out of the possibility of present day comput- 
ers. Moreover, the value of P q at 10 % accuracy level 
can be easily increased by orders of magnitude since ac- 
cording to ([|) t q grows proportionally to y/N. Instead, 
for classical computations, the number of required or- 
bits M will increase exponentially (P ~ 1/M). This 
demonstrates that quantum computations of exponential 
relaxation processes are much more efficient than classi- 
cal simulations with large number of orbits. The pos- 
sibility to efficiently compute the quantum probability 
P q {t) up to very long times t q >> 1/A, allows to numer- 
ically estimate the measure of the integrable component 
fi in the phase space of the classical system. Indeed, for 
/it > 0, the existence of integrable islands leads to an 
asymptotic power law decay of correlations P(t) oc t~ ' 5 
p3[ . Since in our numerical data the quantum proba- 
bility decays exponentially up to P q (t q ) « 5 • 10~ 12 we 
assume that the measure of the integrable component is 
|U < P q (t q ) w 5 • 10~ 12 x (1 ± 0.1) being much smaller 
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FIG. 1. Classical and quantum probability decay for K = 7 
and N/k — 4. The full line shows the fit to the classical 
decay. Dotted lines give the quantum probability P q {i) for 
N = 5001, 20001, 130001 (upper, middle and lower curves re- 
spectively). The lower insert shows the classical probability 
actually computed from M = 9 • 10 9 orbits and the fit is 
shown by the dotted line. The upper insert shows the classi- 
cal (full line) and the quantum (dotted line) asymptotic decay 
for N = 5001. 
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FIG. 2. Mesoscopic fluctuations of the quantum relaxation 
time t q for different system sizes N. The insert shows the 
statistical distribution of fluctuations f which is close to a log 
normal distribution of width 0.22. 
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FIG. 3. The average ratio R(t) = (In ( P q (t) / P(t))} as a 
function of time for different N from iV = 2500 (left curve) 
to TV = 129500 (right curve). The horizontal full line cor- 
responds to P q (t) = P(t). The probability P(t) is given by 
numerical data obtained with M = 9 ■ 10 9 orbits for t < 70 
and by the fit from Fig.l (see text) for t > 70. The left insert 
shows the ratio of the numerically computed classical prob- 
ability P(t) to the fit function. The deviations from the fit, 
for t > 70, are due to statistical errors related to finite M. 
The right insert demonstrates the scaling behaviour of P q (t) 
on the variable t/t q where the t q values are determined by 
condition R(t q ) — 0.1. 
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FIG. 4. Dependence of the quantum relaxation time t q on 
the system size TV = gt c in logarithmic scale. Points refers 
to 10 % deviation level (R = 0.1) while circles refer to 5% 
deviation (R — 0.05). The two dotted lines give the fit 
tq = 1.9iV a41 and t q = 1.5JV 41 respectively. The full line 
gives the theoretical prediction with a = 1/2 (|^). The insert 
shows the 10 % data in a semilog scale. 
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